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Electrons  confined  to  a  double  quantum  well  and  driven  by  a  strong  laser  have 
interesting  properties1'10.  An  electron  that  has  been,  at  a  given  time,  localized  in  one 
of  the  wells  will  oscillate  back  and  forth  between  them.  Grossmann,  Dittrich,  Jung, 
and  H&nggi1'3  have  shown  that  a  laser  with  appropriate  parameters  can  prevent  this 
oscillation  and  keep  the  electron  in  the  well  in  which  it  was  initially  located.  The 
points  {co,E}  (co  and  E  are  the  laser  frequency  and  electric  field  amplitude;  for  which 
this  electron  localization  takes  place  were  called7  points  of  accidental  degeneracy 
(AD).  They  are  given  by  a  simple  equation,  relating  E  to  co,  whose  validity  has  been 
established  under  a  variety  of  conditions1'6-9.  The  Fourier  transform  p(Q)  of  the 
induced  dipole,  which  determines  light  emission  or  absorption  by  the  structure,  also 
has  unusual  features.  If  (co,E)  are  AD  points  |i(Q)  has7'8  intense  lines  at  the 
frequencies  2nco,  where  n  is  an  integer.  This  even-harmonic  generation  is  intriguing: 
the  Hamiltonian  of  the  structure  is  symmetric  and  such  emission  is  forbidden  in  all 
orders  of  a  perturbation  theory  expansion  in  the  field  strength.  Approximate 
equations  for  the  intensity  of  these  harmonics  were  derived  by  Ivanov,  Corkum  and 
Dieterich9'10.  If  the  laser  parameters  are  near  an  AD  point  p(D)  has  a  very  intense 
low  frequency  component7'8  whose  frequency  A  can  be  tuned  by  changing  E  or  co,  or  by 
applying  a  voltage10  across  the  structure.  This  low  frequency  generation  (LFG)  is 
caused  by  the  coherent  oscillation  of  a  major  part  of  the  electron  density  from  one  well 
to  another.  Equations  for  A  are  available9. 

In  this  letter  we  study  the  response  of  an  electron  in  a  double  quantum  well  * 

driven  by  two  lasers  having  the  frequencies  and  co2.  We  show  that  the  two  lasers 
are  capable  of  maintaining  electron  localization  and  generating  low  frequency 
radiation. 

Numerical  solutions  of  the  time  dependent  Schrddinger  equation  for  an 
electron  in  a  double  quantum  well  show  that,  under  the  conditions  used  here,  the 
behaviour  of  the  electron  is  adequately  represented  by  a  two  level  system.  Because  of 
this  we  have  confined  this  study  to  a  two  level  system  whose  parameters  are  chosen  to  » 

represent  a  double  quantum  well.  Analytical  results  are  obtained  by  using  an 
asymptotic  method  developed  by  Dakhnovskii  et  al 9.  They  are  valid  when  v  s  (o)2  - 
coI)/o)2  and  26/hcoj  (2  e  is  the  energy  difference  between  the  levels)  are  small. 
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We  use  the  Hamiltonian: 

H=e{ 1 1X11-12X21)}  -  {1 1X2 1-I2XH ^EUMcvp^Edla*.  (1) 

Here  1 2  )  and  1 1 )  are  the  ground  and  first  excited  states  of  the  electron  in  the  double 
well,  respectively,  and  ox  and  az  are  Pauli  matrices.  The  zero  of  energy  is  halfway 
between  the  energy  levels  e*  and  e2,  and  2e  =  -  £2-  The  electric  field  in  (1)  is  given  by 


E(t)  =  Ej  cosCo^t)  +  E2  coslo^t) 
and  the  induced  dipole  is: 

|i(t)  =  \ii2  <v,t!oxly,t>, 

where  I  vy,t)  is  the  wave  function  of  the  system  at  time  t. 
The  equations  of  motion  for  p(t)  is9 


(2) 


(3) 


dji(x)/dx  =  -  (2e/hco 


4  d^(ti 


)  cos[F(x)  -  F(x*)], 


with 


(4) 


(5) 


F(x)  =  (2  ji^E^o)!)  sin  (t)  +  (2  sin  [(l+v>t]. 

We  use  the  variable  x=031t. 

To  obtain  information  about  ji(t),  for  the  initial  condition  p(x=0)=p12,  we  expand 
the  right  hand  side  of  Eq.  (4)  in  a  power  series  in  (2E/hci)j),  identify  the  largest  terms  in 
this  expansion  and  resume  them.  This  procedure  was  used  previously9’10  for  the  case 
when  the  system  was  driven  by  one  laser.  The  result,  for  the  two  laser  case,  is 

M-(t)  =  p12cos  [  fit)  ] ,  (6) 

with 


fix)  s  (2£/hco1)  I  dxi  Jc[Vei2+  Q22  +  2  eie2COs(v  Tj)  ] , 
JO 


(7) 


and  ej=2  p.12Ej/ho>j.  This  expression  for  y.  contains  only  the  terms  that  were  resumed;  the 
others  are  small,  if  (2e/hcoi)2  «1,  and  are  being  ignored.  The  procedure  gives  the  low 
frequency  part  of  y(t),  hence  the  Eqs.  (6)  and  (7)  can  be  used  to  study  localization  (zero 
frequency)  and  low  frequency  generation. 

If  v  =  0  the  system  is  driven  by  one  laser  having  the  electric  field  intensity 
£=£^£2  and  the  frequency  co  =  oi1  =co2.  For  this  this  case  the  Eqs.  (6)  and  (7)  become 


» 
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|i(t)  =  p.12  cos[At)  (8) 

with 

A  =  (2£/hoo1)J0[e]  (9) 

and  e=2  n12E/hco.  If  J0(e)  =  0  the  leading  part  of  p(t),  given  by  Eq.  (8),  becomes  time 
independent  and  equal  to  p. This  means  that  due  to  the  laser  action  the  system  is,  at  all 
times,  in  one  of  the  states  2* 1/2  [  1 1)  ±  1 2)];  in  either  state,  the  electron  is  localized  in  one 
of  the  wells.  Thus,  the  AD  points  are  given  by  J0(e)  =  0. 

The  induced  dipole  given  by  Eq.  (8)  oscillates  with  the  frequency  A.  By  choosing 
{co,E}  arbitrarily  close  to  an  AD  point  we  can  make  A  arbitrarily  small.  This  is  the  low 
frequency  generation  (LFG)  mentioned  earlier  in  this  article.  LFG  takes  place  because 
we  work  with  laser  parameters  which  are  almost  able  to  localize  the  electron.  As  the 
localization  fails  the  whole  electron  density  drifts  slowly  (on  the  time  scale  A*1)  from  one 
well  to  another;  this  oscillation,  of  practically  the  whole  charge  over  the  large  distance 
between  the  wells,  has  a  very  high  dipole. 

We  examine  next  what  happens  when  the  electron  is  driven  by  two  lasers.  We 
consider  that  the  electron  is  localized  in  one  of  the  wells  if 

p(t)  >0.95  for  x  <  500.  GO) 

Thus  we  require  that  the  induced  dipole  is  very  close  to  the  maximum  value  p12,  which 
can  be  reached  only  if  the  greatest  fraction  of  the  charge  density  is  localized  in  one  of  the 
wells.  The  definition  used  in  prior  work  required  p(x)  =  pl2  at  all  times.  Complete 
localization  for  an  infinite  time  is  possible  only  in  idealized  models  which  ignore  various, 
practically  unavoidable,  incoherent  processes.  For  this  reason  we  prefer  the  definition 
provided  by  Eq.  (10). 

p(x)  depends  on  the  parameters  e^  e2,v  and  2e/hco1.  We  will  explore  this  large 
parameter  space  selectively.  In  Fig.  1  the  dark  regions  show  the  parameter  values,  in 
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the  {v,  e1)  plane,  for  which  the  localization  condition  (10)  is  satisfied.  The  results  for 
negative  values  of  v  are  symmetrical  to  those  for  positive  v.  Since  v  =  0  means  that  only 
one  laser  acts  on  the  system,  the  points  along  the  v  =  0  line  are  given  by  the  condition9 
cos[(2e/hoi)J0[e]  t]<  0.95  for  x<500.  This  gives  small  regions  centered  about  the  roots  of 
the  equation  J0[e]=0,  that  is,  around  the  AD  points  of  the  one  laser  case.  The  size  of  these 
regions  depends  not  only  on  the  physical  parameters  of  the  system  but  on  the  numbers 
used  in  definition  (10)  as  well.  If  we  change  0.95  to  a  larger  number  (but  less  than  one)  or 
500  to  a  longer  time,  the  area  will  shrink. 

The  fact  that  for  v  near  zero  the  electron  localization  takes  place  for  parameters 
located  in  small  regions,  centered  around  the  AD  points  for  the  system  illuminated  by 
one  laser,  is  not  surprising.  However,  the  “fingers”  appearing  in  the  upper  part  of  the 
figure,  for  larger  values  of  v  are  unexpected.  These  regions  are  not  continuously 
connected  to  the  circles  on  the  line  v  =  0.  This  means  that  they  could  not  be  obtained  by 
perturbation  theory  in  v. 

In  Fig.  2  we  show  the  localization  regions  in  the  plane  {e^  e2l.  These  form 
strips  lined  up  along  the  zeros  of  the  Bessel  function  JQ.  The  missing  points  in  the  lower 
left  part  of  the  graph  are  real:  we  have  searched  for  localization  points  there  and  found 
none. 

Finally,  in  Fig.  3  we  show  the  low  frequency  part  of 

dt  e'ifit  W(t-t')  p(t)  |  (11) 

W(t-t’)  is  a  Gaussian  window  function  and  |i(Q)  has  peaks  of  Gaussian  shape  (i.e.  the 
transform  of  the  window  function)  centered  at  the  frequencies  of  the  Fourier  components 
of  p(t).  In  Fig.  3  we  give  only  the  peak  heights  and  positions.  The  parameters  for  which 
the  spectrum  has  been  calculated  correspond  to  the  cross  in  Fig.  1,  which  is  a  point  near 
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the  “finger”  corresponding  to  the  localization  points. 

In  the  system  driven  by  one  laser  the  choice  of  parameters  near  the  AD  points 
lead  to  low  frequency  generation7'9.  This  is  also  true  in  the  two  laser  case:  low  frequency 
generation  occurs  for  all  points  near  those  leading  to  localization.  The  physical  reason 
for  this  is  simple.  At  these  points  the  lasers  are  almost  capable  of  localizing  the  electron. 
As  they  fail  the  electron  density  drifts  slowly  from  one  well  to  another.  This  slow  drift 
corresponds  to  a  large  dipole  oscillating  with  a  very  low  frequency  5.  In  this  particular 
calculation  6  =  1  cm'1.  The  figure  also  shows  what  we  have  called,  in  our  previous 
work7'9  shifted  harmonics,  which  are  Fourier  components  at  the  frequencies  nCo^-cop  ± 
8.  As  the  parameters  approach  the  regions  where  localization  occurs  8  becomes  smaller 
and  ultimately  goes  to  zero;  furthermore,  the  shifted  harmonics  get  closer  and  closer  and 
ultimately  collapse.  We  have  plotted  here  only  the  low  frequency  part  of  the  Fourier 
transform.  At  higher  frequencies  shifted  harmonics  centered  around  n  (al  and  n  <o2  are 
also  present. 

The  low  frequency  part  of  the  spectrum  can  be  understood  from  Eqs.(6)  and  (7). 
Indeed,  fix)  given  by  Eq.  (7)  can  be  expanded  as^13^ 

oo 

fix)  5  (2£/hco1XJo(e1)  J0(e2)x+  £  an  sin(nvx)  /  nv).  (12) 

n=l 

with  a^  =2  (-l)n  J^e^  Jn(e2).  The  leading  term  in  Eq.  (12)  is 

8x  =  (2e/hco1)J0(ei)  Jo(e2)x,  (13) 

which  grows  (when  Jote^)  Jo(e2)  *  0}  indefinitely  with  time,  while  the  other  terms  are 
bound.  Moreover,  the  amplitudes  Jn(ei)Jn(e2)  /  nv  ar€  small  if  v  is  not  too  close  to  zero. 
Introducing  this  in  Eq.  (6),  using  cos(a+b)=  cos(a)  cos(b)  -sin(a)  sin(b)  repeatedly  and 
retaining  the  largest  term  (  the  one  containing  cosines  only)  we  obtain 
cos[fTx)]=cos  [8x]  cos[a1  sin(vx)  /v]  ...cos  [an  sin(nvx)  /  nv]  ... 


(14) 


V 
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Using  in  this  and  the  formula^13^ 

OP 

cos[ansin(nvx)  /  nv]  =  J0(an/nv]  +  2  £  J2m[an/nv]  cos(mnvt).  (15) 

m=l 

we  can  obtain  the  Fourier  components  of  jj.(x>  =  p12  cosfftx)].  We  show  here  a  few  terms 
only: 

^(Typ^  =  ft  J0(an)  cosCSt]  ^0^^  cos[(5+v)x]  +cos[(S-v)x]}/2  +  etc  (16) 

This  equation  is  not  exact.  It  is  based  on  Eqs.  (6)  and  (7),  which  are  obtained  by  collecting 
the  largest  low  frequency  terms  in  )t(x).  For  this  reason  the  formula  does  not  contain  the 
high  frequency  components  of  p(x),  such  as,  for  example,  the  shifted  harmonics  at  the 
frequencies  n  coj.  The  formula  predicts  that  the  amplitudes  of  the  terms  with  frequencies 
(5+v)x]  and  (8-v)x  are  equal.  The  numerical  results  shown  in  Fig.  3  give  nearly  equal 
intensities  as  predicted. 

Taking  the  limit  v  ->  0  in  Eq.  (12),  which  corresponds  to  one  laser  with  the 
intensity  parameter  e=ej+e2  leads  to 

oo 

fix)  ->  (2e/fi<o1){J‘0(e1)  Jote^  £  2  (-l)n  Jn(ei)  Jn(e2) }  x.  (17) 

n=l 

The  factor  multiplying  x  must  equal  A  given  by  Eq.  (9).  Indeed  this  is  the  case  since  the 
expression  in  the  curly  brackets  is  equal  to  J0(e1+e2)14. 

In  this  article  we  have  shown  that  an  electron  confined  to  a  double  quantum  well 
and  driven  by  two  lasers  whose  frequencies  are  closed  to  each  other,  can  be  localized  or 
forced  to  emit  or  absorb  low  frequency  radiation.  In  the  case  when  the  two  frequencies 
are  almost  equal  (i.e.  when  v  =  (o)2~^iVwi  is  almost  zero),  the  present  results  are  a 
simple  extension  of  those  obtained  in  the  case  when  the  electron  is  driven  by  one  laser. 
As  v  is  increased  the  localization  is  no  longer  possible  until  higher  values  of  v  are 
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reached;  then  the  localization  is  again  possible  in  extended  regions  of  the  parameter 
space.  These  regions  are  not  an  analytic  extension  of  those  obtained  for  a  one  laser 
system.  Our  conclusions  follow  from  analytical  results  which  give  the  leading  low 
frequency  contribution  to  the  induced  dipole.  Extensive  comparisons  with  numerically 
exact  calculations  show  that  the  analytical  results  are  accurate  as  long  as  2e/Tuoi  «  1.  In 
real  system  various  incoherent  processes  and  electron  escape  form  the  well  limit  the 
time  over  which  localization  can  be  observed. 

Acknowledgement:  This  work  was  supported  by  NSF  CHE91-12926. 
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Figure  Captions 

Fig.  1.  The  dark  areas  show  those  regions  [given  by  Eqs.  (6)  and  (7)]  in  the  parameter 
space  {v,  e1},  for  which  electron  localization  occurs.  We  used  ej=e2  and  26/hco!  =  0.1  and 
o)l  =100  cm-1. 

Fig.  2.  The  dark  areas  show  those  regions  [given  by  Eqs.  (6)  and  (7)]  in  the  parameter 
space  (ej^)  for  which  electron  localization  occurs.  We  used  2&/haii  =  0.1,  v  =  0.1  and 
0^=100  cm'1. 

Fig.  3.  The  low  frequency  part  of  the  Fourier  spectrum  of  the  induced  dipole  [see  Eq.  (11)] 
for  the  point  in  the  parameter  space  marked  by  x  in  Fig.  1  ^=62=4.5,  2e/ho)i  =  0.1,  v  = 
0.07).  The  spectrum  has  peaks  at  nv±S  where  8=1  cm*1. 
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